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Abstract 

We show that the strong self-interaction of the scalar polarization of 
^ ■ a massive graviton can be understood in terms of the propagation of 

an extra ghost-like degree of freedom, thus relating strong coupling to 
the sixth degree of freedom discussed by Boulware and Deser in their 
Hamiltonian analysis of massive gravity. This enables one to understand 
the Vainshtein recovery of solutions of massless gravity as being due to 
the effect of the exchange of this ghost which gets frozen at distances 
larger than the Vainshtein radius. Inside this region, we can trust the 
two-field Lagrangian perturbatively, while at larger distances one can 
use the higher derivative formulation. We also compare massive gravity 
with other models, namely deconstructed theories of gravity, as well as 
DGP model. In the latter case we argue that the Vainshtein recovery 
process is of different nature, not involving a ghost degree of freedom. 
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1 Introduction 



The construction of a consistent theory of massive gravity has proven to be an 
extremely difficult task. The usual Pauli-Fierz approach is to add a quadratic 
mass term to the quadratic Einstein-Hilbert action for a perturbation (the would 
be massive graviton) over a flat background |Jj. At this level, the action for the 
massive graviton is consistent in the sense that it does not propagate ghosts or 
tachyons. Nevertheless, it is incompatible with experiment because it propagates 
an unwanted extra scalar degree of freedom (DOF) that couples to the trace of the 
energy momentum-tensor j5| (further referred to as the vDVZ scalar). It leads to 
physical predictions different from the massless theory, even in the limit where the 
graviton mass goes to zero. This is true for instance for the light bending around 
the sun, and is known as the van Dam-Veltman-Zakharov (vDVZ) discontinuity j2]. 
An other question raised by the quadratic Pauli-Fierz theory is the link between 
the background flat metric and the graviton. This should be answered by a proper 
nonlinear theory of massive gravity. However, attempts to go beyond the quadratic 
Pauli-Fierz action are raising even more questions. For example, Boulware and 
Deser j3j introduced an extraneous metric, so that the theory they considered can 
be thought of as some kind of bimetric theory jl] , the dynamics of one of the two 
metrics being frozen, and the coupling between the two metrics being such that it re- 
produces the quadratic Pauli-Fierz action for small perturbations. In the following, 
we will refer to this type of theory, quite loosely speaking, as massive gravity, note 
however that such a nonlinear completion is not unique. In a particular example of 
such a theory, simply obtained by adding the quadratic Pauli-Fierz mass term to the 
full Einstein-Hilbert action, a standard Hamiltonian treatment shows that the num- 
ber of propagating DOFs is six (instead of the five DOFs described by the quadratic 
action), of which the sixth DOF has a ghost character This increment in the 
number of the propagating DOFs with respect to the quadratic theory, as well as 
the unboundedness from below of the Hamiltonian, was shown to persist in a more 
general class of bimetric theories [3] . This could have consequences at the classical 
level already (see e.g. [U|). On the other hand, it was argued in [7J|H| that the vDVZ 
scalar has a strong self interaction on a flat background at an extremely low scale 
A = (Mpm 4 ) 1 / 5 . This, at the classical level already, has dramatic consequences. 
Indeed, Vainshtein j^j first noticed, that nonlinearities in massive gravity dominate 
the linear terms at distance scales smaller than ry = (GM/m 4 ) 1 ^ 5 (the Vainhstein 
region) for a classical source of mass M. He then proposed that these nonlinearities 
cure the vDVZ discontinuity within this region, due to a nonperturbative resumma- 
tion of the source expansion [0] (see also [7] ) . This " Vainshtein resummation" seems 
to be problematic for massive gravity ^U], if one follows the resummation beyond 
the first terms, but there remains the possibility it does work for more sophisticated 
models (see e.g. [7]). Nonlinearities could also cure some of the above mentioned 
pathologies by selecting a non asymptotically flat vacuum [3] . 

In this paper, we will only consider expansions of massive gravity over a flat 
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background. Namely we want to discuss the relation between the strong coupling 
and the, seemingly unrelated, ghost appearing at nonlinear level in massive gravity. 
That there could exist a link between these two issues is suggested in part by the 
following reasoning. Consider e.g. a classical non-relativistic source. There it seems, 
if one follows Vainshtein's reasoning, that within the Vainshtein region, one somehow 
compensates for the extra attraction exerted by the vDVZ scalar, that is responsible 
(in a perturbative sense) for the discontinuity (this compensation being at worst only 
valid far from the singularities discussed in [J0J). One may suspect that there is a 
ghost degree of freedom responsible for this cancellation 5 . A natural candidate for 
this ghost seems to be the "sixth degree of freedom" mentioned above, but it is 
not clear why this ghost should only work within a certain distance scale. We will 
discuss in this paper how one can reformulate the strong coupling as a ghost problem 
in the theory, and our analysis will clarify the above features. 

Although we will focus here on massive gravity (as defined above), we would 
like to motivate this work by more recent developments. Indeed, massive Einstein 
gravity is just one example of a theory that modifies gravity at large distance-scales. 
In recent years, a number of different models have been proposed to modify gravity 
in the infrared. In particular, the Dvali-Gabadadze-Porrati (DGP) braneworld 
model (also known as brane induced gravity), modifies gravity at large distance 
while it can give an alternative way of producing the observed late-time cosmic 
acceleration . DGP gravity shows certain similarities with massive gravity. More 
specifically, it exhibits a vDVZ discontinuity at linearized level, and also has a related 
strong coupling [3 HE1 HBj, the exact consequences of which has been subject to a 
debate [7J HEl H3 HBj- With the link between strong coupling and ghosts in 
massive gravity, one would hence suspect the DGP model also to contain ghosts 
when properly analyzed (i.e. at the nonlinear level). As we will argue, the situation 
in DGP is quantitatively different from massive gravity, in that the leading operator 
that grows strong in the UV has dimension seven, as opposed to massive gravity, 
where this operator has dimension nine. This fact and the particular tensor structure 
of the operator describing the scalar content of the graviton, as found in [To^ ITT] , 
suggest that the issue is more subtle and differs from the case of massive gravity. 

This paper is organized as follows. In the next section we will discuss how, 
in massive gravity, the strong coupling can be reformulated as a ghost problem. 
Namely we will argue that there is a propagating ghost DOF, which appears only 
at the cubic order in the perturbation theory over a flat background, and which 
is responsible for the cancellation, a la Vainshtein, of the attraction exerted by 
the vDVZ scalar around heavy sources. We will show how the ghost formulation 
allows a perturbative treatment of massive gravity within the Vainshtein region (so 
that in our reformulation, the theory stays weakly coupled in this region). This 
perturbative ghost generates the same leading order corrections as found in the 
Vainshtein resummation as discussed in jHHTJ, while it freezes out at large distance, 

5 Such a mechanism was also found in the GRS model [11] . responsible for the instability of the 
model P2]. 
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where one is only left with the vDVZ scalar. This relates the old Hamiltonian 
formulation to the effective field theoretical formulation. In section 3, we want 
to understand how the appearance of ghosts in massive gravity is related to the 
breaking of accidental symmetries, present in the linearized theory, at nonlinear 
level. There we will discuss the link between the ghost discussed in section 2, and 
the "canonical" sixth ghost DOF as discussed by Boulware and Deser. This section 
will further give a "geometrical" meaning to the ghost problem, which we further 
explore by considering discretized gravity in section 4. There, we also compare 
massive gravity with the DGP model. We suggest that the situation in the latter is 
very different as far as the relation between strong coupling and ghosts is concerned, 
namely because the operator that grows large at the lowest energy scale in the DGP 
model, is not associated with propagating ghost modes. We conclude with some 
final remarks on our work and other IR modifications of gravity. Note that our 
discussion will be purely classical, and we will not be concerned with the issue of 
the quantum consistency of the theories considered. 



As recalled above, the quadratic mass term for a spin two excitation is uniquely 
defined by demanding the absence of ghost and tachyonic modes, and takes the 
Pauli-Fierz form pQ. Beyond quadratic level, completion is not uniquely defined. 
Two possible choices are for example given by the following actions 



where h^ u = g^ u — g^J, g\J being some extra metric field 6 . The above actions 
define bigravity theories, where the dynamics of g$ has been frozen 7 . Note that the 
above theory is invariant under 4D diffeomorphisms (in the same sense as bigravity 
theories considered in [3] are) and arise naturally in the process of "deconstructing" 
gravity jKJEi] . The action (0) was used by Boulware and Deser in their Hamiltonian 
analysis of massive gravity [3], while action (J2J) was used by Arkani-Hamed et al. in 
the ref. jS]. In this last paper, an action for the scalar polarization of the massive 
graviton was obtained by introducing the "hopping" fields Y a (x), and replacing 

6 and g^ is the inverse of gffl 

7 This theory is in the "Pauli-Fierz universality class" , to use the phrasing of ref. |S] . 



2 Ghost or Strong Coupling? 
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by given by 

H iaf = g iaf -g^Y)d l y a d v Y^ (3) 

Taking gj$ and g^ v to be flat, and expanding Y a as Y a = x a + ir a , one obtains 

Hfj, u = + tt^ u + 7T U)fl + 7r a)M 7r". (4) 

Focusing on the scalar mode of the Goldstone vector tt^, ir^ = <9 M 0, one gets (after 
an integration by part) the cubic action for 

2M 2 P m 2 J d 4 x ((D0) 3 - {U < t>){& 1 d v (j>){d tl d v (j>)) . (5) 

Note that, due to general covariance of the Einstein- Hilbert action and to the partic- 
ular tensorial structure of the Pauli-Fierz mass term, there is no kinetic term arising 
directly for from the above procedure. Rather, the Goldstone scalar obtains a 
kinetic term only through mixing with h^. This gives a m 2 dependence (and pos- 
sibly very small coefficient) to the scalar kinetic term ~ Mj,m 4 0[Z]0. After a proper 
diagonalization jH], one obtains the following action for the canonically normalized 
scalar (we will still denote 0) 

C = ^ + ^(D0) 3 - ^(□0)(^0)(d A A0) - ^0T, (6) 

where T is the trace of the energy-momentum tensor, and the energy scale A is given 
by 

A=(m 4 M P ) 1/5 . (7) 

This Lagrangian contains strong cubic interaction terms for the Goldstone field, with 
dimension nine operators growing strong at the scale A. One can check that the 
cubic terms are the dominant terms in the effective Lagrangian, in the sense that 
they grow strong at a much smaller energy scale than any other interaction terms 
[8J, which makes this Lagrangian a useful description of the low-energy dynamics. 
This amounts to taking the limit Mp —>■ oo, T — >• oo and m — > with A = est 
and T /M p = est in the original Lagrangian (J2J). The limiting procedure eliminates 
all other self-interaction terms of as well as mixing terms of the Goldstone scalar 
with other spin components. A similar limit has been considered in in the case 
of DGP gravity. Note also that, had we applied the same procedure from action (JTJ), 
we would have obtained the same action as in © with a global minus sign in front 
of the cubic derivative interaction terms. This means that the following discussion 
will apply for both cases, and indeed, as we will argue, much more generally. 

To investigate the properties of the 0-sector, and with no fundamental differences 
in the conclusion, we simplify our discussion here, by omitting the third term in the 
Lagrangian © (in the Appendix we explain how to extend the treatment below for 
the full Lagrangian). So the starting point in our discussion will be the following 
Lagrangian 
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This Lagrangian obviously implies an equation of motion of fourth order, reading 

□0 + ^n((D0) 2 )-^ = o. (9) 

This means, if one thinks in terms of the Cauchy problem, that one can expect 
this action to describe two, rather than one, scalar DOF. As is well known, similar 
conclusions are reached in higher derivative scalar field theories (see e.g. [20 J or 
higher derivative theories of gravity (see e.g. [21] )• One can typically reformulate 
a higher derivative theory of some fundamental fields as a standard two-derivative 
Lagrangian of the fundamental fields plus extra degrees of freedom, with certain 
interaction terms, and which encode the extra derivatives appearing in the initial 
action. It is also known that some of these extra DOF, quite generically, have a 
ghost character (see e.g. [22]), so that one should be able to reinterpret, classically, 
the strong coupling discussed in [7J [S] in terms of an extra propagating ghost DOF. 
This is what we do in the following. Note that one could question the consequence 
of this because we consider here a truncated theory. For example, if one looks at a 
truncated derivative expansion of some perfectly safe underlying theory, one could 
conclude erroneously that the latter is sick. Here, however, the situation is quite 
different because of the fact the Lagrangian with started from has some well defined 
range of applicability, as was recalled previously. 

We thus introduce a new field A, modifying the Lagrangian (JHJ) into 

C eq = + ^(D0) 3 + F(\, D0) - ^0T. (10) 

The equations of motion for A and <fi are given by 

□0 + ^□((□0) 2 ) + UF<M - = 0, (11) 

and 

F (1 ' 0) = 0, (12) 

where F™' means the derivative of F, i times with respect to its first variable and 
j times with respect to its last one. We then ask those equations of motion not to 
contain derivatives of order higher than two and to be equivalent to equation Q. 
A suitable ansatz for F is given by 

F(A, D0) = ^A 5 / 2 A 3 + A 2 D0 - ±r (D0) 3 , (13) 



leading to the equation of motion 



□0 + DA 2 -— T = 0. (14) 

M P v ' 

/A 5 / 2 \ 

a(— -A + □</>) =0. (15) 
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Those, for a non- vanishing A are equivalent to the equation of motion Q. Note that 
the phase space of solutions to the equations (|14H15|) is in fact larger than that of 
solutions to Q, because the former also includes solutions to the equation (obtained 
for a vanishing A) □</> + T/Mp = 0. It means that the equivalence we are talking 
about here has to be understood as an equivalence restricted to a suitable subset of 
solutions. Similar restrictions also occur, e.g., when showing the equivalence between 
a f(R) theory of gravity (/ being some arbitrary function) and a scalar tensor theory, 
where the equivalence only holds between solutions defined away from the critical 
points where the second derivative of / vanishes [2T] (see also (2H1)- The choice (fT3Jl 
leads to a Lagrangian where the field A 2 does get a kinetic term through a mixing 
with <p. Defining then ip as <p = tp — A 2 , one obtains the Lagrangian 

C eq = -cpDcp - -X 2 D\ 2 + — ^A 3 A 5/2 - — ^pT + — A 2 T. (16) 
q T y 2 3^ M P y M P v ; 

This can be then rewritten as a quadratic action defining ip as ip = A 2 , so that C eq 
becomes 

£e „ = \ vap - 1** - ^ W - -&r + -fcr, (17) 

where e = ±1 is the sign of □(<£> — ip). Notice above the non-analytic form of the 
potential for ip, as well as the fact that ip is a positive definite field variable. The two 
different possible choices for e in the above equation indicates that the equivalence 
between the dynamics obtained for <fi from the Lagrangian (jHJ) and that obtained 
from the Lagrangian (JT7J), with a specified value of e is restricted to a set of solutions 
to the equations of motion having the corresponding sign of □</> = — ip)- If one 
consider e.g. the Cauchy problem associated with the equation of motion 0, and 
if one wants to derive the solution to this Cauchy problem with the equivalent 
Lagrangian (f*T7j). one should choose the sign e corresponding to the initial Cauchy 
data provided (of course the sign of □</> can change along the initial surface so that 
the equivalence we are talking about here is really a local property). In any case 
one can also use the Lagrangian (fTB^l where this sign issue does not arise, but as we 
saw, this also leads to solutions of the equations of motion not contained in those 
of ©. 

Thus, we rewrote a Lagrangian of a field with nonrenormalizable interactions as 
a Lagrangian of two degrees of freedom; one free field, ip, corresponding to the vDVZ 
scalar at large distance from a source, and one ghost ip, with a relevant interaction 
term. 

Let us for now consider a point like source. Then, within a certain distance 
scale to the source, the ghost exactly cancels the vDVZ field, up until when the 
ghost freezes out due to its self-interaction, leaving one propagating DOF. For our 
discussion, it is convenient to consider the dimensionless fields (p = (p/Mp, and ip = 
ipjMp which correspond with no further normalization factor to the dimensionless 
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massive graviton h^ v . The equations of motion for those dimensionless fields read 

Hip = GT (18) 

D ^ + e ^l/2 = GT; (19) 

where G = j^t. Considering a classical source with T = — M5 3 (x), one finds that ip 
is given by the usual Newtonian potential (with an appropriate sign) 

= ^, (20) 
r 

where r is the distance to the source. If one then considers the expansion of the 
solution for ip around ip^ = <p(°\ as ip = ip^ + ip^ + one finds that ip^ obeys 
to the equation 

□^(D +e ^!^(o)y /2 = o, (21) 

which is solved by 

4 



^« = -e -m 2 VGMr 3/2 . (22) 

15V3 

This matches the correction 8 obtained in jS] (see jjj). Here we traced the origin 
of these corrections back to the self-interaction of the perturbative ghost ip, which 
cancels the contribution from the vDVZ scalar ip at small distances from the source. 
Notice that ip^ becomes of the same order as ip^ at the distance ry = (GM/m 4 ) 1 ' 5 , 
which is the Vainstein radius j^j. Going back to the original field = ip — ip, we 
thus see that, inside the Vainshtein region <p^ = and 



{1) = e =m 2 VGMr 3/2 r < ry. (23) 

15v3 



Outside the Vainshtein region, one cannot trust perturbation theory for the ghost, 
but in this region, one can simply use the original Lagrangian © for the scalar field 
perturbatively, since the self-interaction of the scalar is small in this region. The 
equations of motion for this Lagrangian read: 

□0 + |rD(D0) 2 - ±n (,9^0(9^0) - ((D0)^^0) - T/M P = (24) 

If we now look for an expansion of the form = ( - ^ + < - 1 - ) + ... with < - 1 - ) <C 0^°\ 
we see that cp^ obeys 

□0(°) = T/Mp (25) 

so that 0*- ** is given by 

i r 

(26) 



rMr 



8 One can verify that Dip obtained from this solution has indeed the same sign as e. One should 
choose the sign of e so that the solution has the right asymptotic behavior. 
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for a point like source. is now given by 

□0« = J_D (d^^d^W) (27) 

(the terms in d^ ) vanish outside of the source...) This is solved by 

M 2 1 

(1) oc— ^— — r>r>, (28) 

which is the form of the correction expected (recall that the dimensionless metric 
fluctuation are given by <f>/Mp). 

To summarize, we rewrote the Lagrangian for massive gravity in which the vDVZ 
field gets strongly coupled at a scale A, as a system of two fields that do not have 
strong coupling, but of whom one is a ghost. In this ghost formulation, we can trust 
the Lagrangian perturbatively within the Vainhstein (or strong coupling) region, 
while outside of the Vainshtein region, we can use the original higher derivative 
interaction perturbatively. 



3 Symmetries and Hamiltonian analysis of Mas- 
sive Gravity 

The " Goldstone" formalism that we used as a starting point in the previous section is 
especially useful to focus on the sick scalar sector of massive gravity. Here we relate 
our discussion to the original analysis of Boulware and Deser, who first discussed, in 
the "unitary gauge", the unboundedness from below of the Hamiltonian of massive 
gravity [3]. For that discussion, our starting point will be the action (0), even if, as 
we already stressed, the Boulware and Deser findings have been shown to be valid 
for generic bimetric gravity theories (and thus massive gravities as considered here) 
by D amour and Kogan. 

Let us first recall what is going on, from the Hamiltonian point of view, if one 
only retains the quadratic part of the action (0) for the "graviton" h^ v . That is 
to say, one sets in (JTj) the background metric g$ to be the Minkowski metric, and 
expands the Einstein-Hilbert action (first line of (JTJ)) to quadratic order in h^ v . One 
obtains then the Pauli-Fierz action, which is ghost-free and propagates five degrees 
of freedom. The counting of DOFs goes as follows. One start from a symmetric 
tensor, the graviton, which has 10 components. However it turns out that neither 
/loo nor h Qi (with i = 1,2,3 being spatial indices) are dynamical degrees of freedom. 
This is a consequence of the same properties of massless gravity which shares the 
same kinetic terms with the theory considered here 9 . Namely /loo an d h 0i appear 

9 For the same reason, the A component of a massive Proca fields does not propagates as is 
the case for the massless photon 
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as Lagrange multipliers in the kinetic term obtained from expanding the Einstein- 
Hilbert action. This means that out of the 10 initial DOF one is left with 10-4= 6 
DOF at this stage. Let us then consider the mass term, here h 00 and h 0i play very 
different roles. Indeed the mass term reads 



M P m J d x (hijhij - 2h 0i h 0i - huhjj + 2huh m ) (29) 

where it appears that /i o is a true Lagrange multiplier for the theory considered, 
since it enters linearly both in the kinetic part and mass term of the action. As a 
consequence the hoo equation of motion generates a constraint reading 

(V 2 - m 2 )/^ - h m = 0, (30) 

which enables to eliminate one more DOF leaving 5 propagating DOF in the Pauli- 
Fierz action. Note that the equations of motion for the h 0i , which enter quadratically 
in the mass term, do not eliminate other degrees of freedom (in contrast to the 
massless case), but rather determines those variables in terms of the others. 

We now turn to fully nonlinear gravity. In the massless case, that is to say Ein- 
stein General Relativity, formulated in the Hamiltonian language (ADM formalism 
[21]), the lapse and the shift field N and shift N l fields are Lagrange multipliers 
associated with the reparametrization symmetry of the Einstein-Hilbert action 10 . 
The latter are defined as 

iv=H/T 1/2 (si) 

and 

Ni = goi , (32) 

in terms of the component of the metric g^. They generates constraints related re- 
spectively to the time and space reparametrization symmetries of the action. Those 
constraints eliminate 4 degrees of freedom, out of the 6 remaining (10-4) leaving the 
two propagating polarizations of the massless graviton. The addition of a "mass 
term" such as the one of the second line of ([!} changes however dramatically the 
character of iV and N l . Indeed, the action (0) reads in the first order formalism 



M 2 P 



J d 4 x { (n ij gij - NR° - JVii?) 

(hijhij - 2NiNi - huhjj + 2ha (l-N 2 + N k g kl Ni) ) } , (33) 



-m 2 



where 7T 4 - 7 are the conjugate momentum to and R° and R l are respectively the 
Hamiltonian and momentum constraints of massless gravity (generated by the lapse 
and shift fields). The dramatic observation made in ref. [3] is that now, neither iVj 
nor iV are true Lagrange multiplier of the full non linear massive gravity. Thus, the 
number if propagating DOF is now 6 rather than five. It is also remarked in ref. [3] 

10 that is to say, the Einstein-Hilbert action, in its first order form, is already in the 
"parametrized" form, see e.g. |24l ES] 
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that the full reduced Hamiltonian is unbounded from below as we now recall. For 
this purpose, we introduce the variable n defined as 



N = 1 + n (34) 

and rewrites the action (|33|) as 

, J d 4 x{ (ir ij gij - (1 + n)R° - N t IV) 

(hijhij - 2NiNi - h,;h u - 2h u (2n + (3n 2 - aN k g kl N t ) ) } , (35) 



M 2 P 



-m 



with (3 and a equal to one. For future reference, we will however keep explicitly 
the (3 and a dependence upon the process of reducing the Hamiltonian. Notice in 
particular, that upon the substitution iVj — > h Qi and 2n — > —hoo, the mass term ()29|) 
is obtained from the above expression by taking (3 and a to vanish. The iVj and n 
equation of motion read respectively 

R = 4m 2 (r] ij — ah kk g ij ) Nj (36) 
R° = Am 2 hii(l +(3n). (37) 

Those can be used to extract n and N l 

Nj = -^(ri-aJwtflt (38) 

n = A ol 2 ( R ° ~ 4m Xi) » ( 39 ) 

Inserting those expression in action (J33J), we get the reduced Lagrangian (in first 
order form), 



M 2 P 



J d A x ^iT^gij - m 2 (hijhij - hahjj) - ^2 Rl {v ~ ah a3)im Rm 

1 (it! ) 2 -^4" ( 4 °) 



8m 2 f3hii v ' 13 

One reads immediately from the above expression the Hamiltonian of the system, 
and discovers that it can have arbitrary sign and absolute value, since it is true in 
particular for the first term of the second line of ([40)1 . This expression (J4*U|) shows 
in a very clear way that the number of propagating DOF is now six. Note however 
that the Hamiltonian is singular when one makes f3 to vanish. This manifests the 
fact the sixth DOF only arises from the integration of n, as well as that we have 
chosen to integrate the constraint ()39|) by extracting n, which is not possible if (3 
and ha vanish. So, in this language, the decoupling of the sixth DOF is subtle and 
the above Lagrangian can be considered analogous to the Lagrangian (|17|) where 
the two degrees of freedom appear explicitly. 
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We would like now to relate the unboundedness of the Hamiltonian for the extra 
propagating DOF to the analysis done in the section |2] of this paper. We will discuss 
this issue in the covariant formulation, following first the discussion of [3]. 

Starting from the action (fl|). one obtains the following equations of motion for 
the metric g^. 

g^{g = r, + h)- 2m*^Mh af} (gftg* - gfigfy = GT^ (41) 

where Q^ v is the Einstein tensor for the metric g. In this formulation we can again 
see how at linearized level only five DOF appear. Indeed, by taking the divergence 
and trace of the linearized form of the equations of motion, one obtains at the linear 
level 

d^hT - h-rT) = 0, (42) 

h = --^T. (43) 
orrr 

The first equation implies that the linearized curvature R L is zero, R L being given 
by 

R L = 28^ {h^ u - d*hl) . (44) 

The second equation states that the trace is not propagating. The latter constraint, 
obtained by taking the trace (with respect to g^ v ) of equation (JHJ), generalizes to 
the full nonlinear case into 

~R ~ 2m 2 ^V^ (<oWo) - ^0)4) = GT ( 45 ) 

which, involving second derivatives, has to be treated as an equation of motion. This 
is indeed where the sixth ghost degree of freedom comes in the covariant formulation. 
Now the first constraint in eq. (j42j) has a nonlinear generalization which is obtained 
from the Bianchi Identities for the metric g. Indeed taking the covariant derivative 
with respect to g^ v of (JHJ, one obtains 



■j=z \9(o)9(q) ~ 9(0)9(0) ) I - U. (4b) 



-9 

Taking the ordinary divergence of this equation allows one to write the linearized 
curvature in terms of higher order terms in h^. Schematically the above equation 
reads 

d v d»[W v -rTh] + nh 2 + 0{ti i ) = 0, (47) 

yielding Ri(h) = Oh 2 up to terms of 0(h 3 ). When plugging back in the trace 
equation (J43j) . we get, again schematically, at quadratic order in the field equations: 

Oh 2 + m 2 h = GT. (48) 

If we replace h 2 — > if), we see that this equation has exactly the same form as the 
equation for the ghost of the previous section eq.(fT§j). 



12 



4 A comparison with other models 



4.1 Deconstructed Gravity 

In the previous section, we saw that the Pauli-Fierz mass term was such that the 
equivalent of the R° = constraint, associated with time reparametrization of Ein- 
stein (massless) General Relativity, was still present in the massive theory at the 
level of the quadratic action and is given by equation (J30|) . Thus, somehow, the time 
reparametrization symmetry is still present in the quadratic Pauli-Fierz action and 
eliminates two degrees of freedom. However it disappears from the fully nonlinear 
massive gravity, as we just reminded. It is interesting to remind a similar property 
arising in the process of deconstructing gravity [HI 1213 CHj, where geometry along 
one (in the simplest case considered here) dimension is given up by discretization. 
The starting point is here the five dimensional Einstein-Hilbert action reading 



J d A xdy^-gM {R + K^Kap (g^g aP - g» a g vP ) } , (49) 



where we used an ADM-split of the 5D metric along the fifth space-like dimension 
to be discretized (see [19]). One then latticizes the fifth dimension, with S lattice 
sites (labelled with latin indices (i), (j), ...) with spacing a (we will assume periodic 
identification of the end point of the lattice for simplicity). We obtain the discretized 
action 

S = Ml 5) aJ2 J d"x^g~^M (i) [R(g®) 

l (50) 
— (A c g (l) )^(Acg {i )) a( 3(g^g^ - 9$ 9$)] , 



(0 



where we took the natural discrete definition of the Lie-derivative 

a 

In this expression W is the transport operator from site i + 1 (located at coordinate 
ya+i) along the fifth dimension), to site i (located at coordinate = yu+i) ~ a 
along the fifth dimension). It is given by the Wilson line for transport along the 
fifth dimension (in analogy with non-abelian gauge theory) 

W(y',y)=Pexp f dzN»d„, (52) 

with iV 7 ^ are the the shift fields of the ADM-split used in this discretization. The 
hopping field Y^ 1 (analogous to those of equation (j3J|) considered in [H] can be ex- 
plicitly constructed out of W as ?/o; = W(y,yo)(x)- 
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It turns out to be convenient to work in the Einstein frame for the metrics on 
the different sites. Namely, we perform a Weyl rescaling g$, = exp ^— ^ j 7^, and 

Mi = exp(0( i )/v / 3), such as the action (joTIj) now reads 

V ^ 2 (53) 

with 



0{2 = i{Ar7®-7{2^^ (54) 



and Mp = M%^a. Eventually, we want to discuss the symmetries and counting of 
DOF of this action at the quadratic level for some fluctuation over a flat background. 
We thus expand the different fields as follows 



* = "W + £e = ^(M + l) = ** + ^„ (55) 

and define furthermore the discrete Fourier transform of the discrete fields as 



.1 

This gives the following quadratic action 



T {k) = £ -^e-*"*/". (56) 



k 

k 

+ V - sin 2 — ( f f ji« - ^<3«) - 2a9{ ^ ) ( ( h*^ - M<**(*A 



.*(fe) 



2ad( a n 



e -i2 

As we will now argue, this action describes, leaving aside zero modes, a tower of 
massive spin-two fields with a mass spectrum given by 

Indeed, if one concentrates on the massive modes, the symmetries of this action are 
as follows. First one has the following Stuckelberg symmetry acting at each mass 
level, namely 

( C i2ivk/N _ -1 \ 

^ = 2d ( ,e } \ H fc) = ^ (59) 

(X 
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so that the vector fields are the Goldstone bosons which get absorbed by the (5 — 
1) gravitons that become massive. This symmetry was fully expected from the 
symmetry of the original action (J5U|) since it comes from its invariance under the 
product of all the 4D diffeomorphism invariance on each site. It eliminates 8(5'— 1) 
DOF out of the 105 + 4(5 - 1) + 5 = 155-4 DOF present in the action and 
corresponding respectively to the 5 4D metrics at each sites, the 5 — 1 hoping fields 
Y in between adjacent sites and the 5 lapse functions on each site. However, there 
is an extra accidental symmetry acting at the quadratic level (in the action), the 
latter reads 

= v»»f {k \ # W = VSf {k \ SUM = l _ e a i2 ^ k/N d,f {k \ k ? 0. (60) 

and eliminates an extra 2(5 — 1) DOF. We see that this leaves 5(5 — 1) propagating 
DOF (not counting zero modes) corresponding to the announced (5 — 1) massive 
gravitons. Interestingly, this symmetry is inherited from the reparametrization in- 
variance along the discretized direction which is however broken by the discretiza- 
tion. This is thus analogous to what was shown to happen for the Pauli-Fierz theory 
with the time reparametrization invariance. In particular it is not expected one can 
extend this symmetry at the nonlinear level so that extra degrees of freedom will 
start propagate there. 

4.2 DGP gravity 

The DGP model is a four dimensionally covariant braneworld model that is closely 
related to massive gravity [T3]. The DGP setup describes a three brane in a five 
dimensional flat bulk, on which a large Einstein-Hilbert term is present. The action 
defining the model is taken to be 

S = M 2 P j y/\jj\ Rd A x + J ^/\i\Rd 5 X, (61) 

with Mp 3> M* and where a usual Gibbons-Hawking term taking care of the brane is 
implicit. What is important for us is that, from the point of view of an effective four 
dimensional brane observer, the DGP action produces infrared modified Einstein 
gravity in a way related to massive gravity. Indeed, in a particular gauge, the 
fluctuations around Minkowski vacuum are described by a Pauli-Fierz alike effective 
action for a massive graviton with "running mass" m 2 = p/L^gp, where p is the 
graviton momentum and L^gp = M 2 /M^ is the DGP length scale. This scale was 
found to be the crossover scale between a small distance four dimensional behaviour 
of the Newtonian potential exchanged by non-relativistic sources on the brane, and 
a large distance five dimensional behaviour |^. As in Pauli-Fierz theory, DGP 
theory propagates gravitons with five polarizations at linearized level. As such, 
one might expect that the model shares the same difficulties as massive gravity, 
especially the strong coupling and ghost problem. Indeed, the former property of 
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DGP was made explicit in and JH] in confirmation of the work j7j discussing 
the appearance of a Vainshtein scale in the model. The exact consequences of this 
strong coupling is subject to a debate [T71[TB|. Moreover exact solutions |2Z] as 
well as approximate ones [28J indicate that a Vainshtein mechanism is at work in 
the DGP model, allowing to recover solutions of Einstein General Relativity (see 
however [2*3 )■ In the light of our previous discussion of massive gravity, one might 
suspect that this Vainshtein "resummation" is, also in the DGP model, due to a 
dynamical ghost DOF at work. However, with the DGP model being 4D covariant, 
one should investigate this issue in detail. There is one obvious first qualitative 
difference in between the two theories. Indeed, the operator in the DGP effective 
action, which grows large at a low scale (the strong coupling scale), is a dimension 
7 operator ^3] , as opposed to the dimension 9 operator found in the the Pauli-Fierz 
action. The effective theory for the longitudinal mode of the graviton in DGP is 
obtained by decoupling in a similar way to what was recalled in section E]for massive 
gravity. It reads [T7j 

1 7rT 
C dgp = 3nDn - — (d,7r) 2 Dn + —. (62) 

1V dgp Z1VI P 

The scale at which the 7r-mode gets strongly coupled is given by A dgp = (M p / L 2 ,^) 1 / 3 . 
One might think that the higher derivative operator appearing here-above also cor- 
responds to extra degrees of freedom in an equivalent formulation. However, the 
equation of motion for this scalar are easily seen to be: 

GDtt - -n^(cVU) 2 + -n^) 2 = (63) 

^dgp Jy dgp 4iW P 

Although this is a nonlinear equation of motion, it is a second order differential equa- 
tion for the mode it, unlike the analogous equation for its Pauli-Fierz counterpart. 
Hence, as seen from a Cauchy problem perspective, the equation describes only 
one propagating degree of freedom. In J7|, the geometrical meaning of this equa- 
tion was clarified by noting that in the full theory it descends from a combination of 
the Gauss-Codazzi equations and 4D Einstein equations on the brane, which, by the 
relation ~ —ji—(d^d u 7i) makes eqn.fpjB*)) an algebraic equation for the extrinsic 

dgp 

curvature of the brane. Considering now a point like non-relativistic source, as done 
at the end of section |21 one finds that inside the Vainshtein region, the solution for 
7T is given at dominant order by the solution to the the quadratic part of equation 
(JB3|), that is to say the equation obtained by dropping the first term in the left hand 
side of (Jfillj) [17J. Thus, in order to estimate the contribution of tt to the Newtonian 
potential around the source, one may say that 7r obeyes in this region an equation 
similar to equation (|2*T|) obeyed by ip^ 1 ', with the role of ip^ being played by T/M p . 
This yields indeed corrections to the Schwarzchild solution that goes (correctly) as 
(r/r*) 3 / 2 , whgj-g r * = (j-^QM) 1 is the Vainshtein radius for the DGP model. How- 
ever, here there is no need for the cancellation of an extra contribution coming from 
another DOF. 
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Note this does not tell us what is the nature of the small fluctuations around some 
background solution. In particular, in ^7j it is shown that the kinetic term for this 
scalar is enhanced significantly in the Vainshtein region, which incidentally means 
that the 7r-mode interacts more and more weakly as one approaches the source. 

5 Conclusion 

In this work, we compared the Goldstone formalism for massive gravity jH] , to the old 
Hamiltonian approach of Boulware and Deser [3j . In the first approach it was found 
that the scalar polarization (the vDVZ scalar) of the graviton acquired a strong cubic 
interaction, while in the second approach it was shown that, at nonlinear level, mas- 
sive gravity propagates a sixth degree of freedom with negative energy. We showed 
that one can reinterpret the strong coupling of the vDVZ scalar as the propagation 
of a ghost, in agreement with the Boulware-Deser finding. The ghost can then be 
seen as responsible for the cancellation of the attraction exerted by the vDVZ scalar 
at distances around a non-relativistic source smaller than the Vainshtein's radius. 
Inside this region, we can trust the two field Lagrangian perturbatively, while at 
larger distances one can use the higher derivative formulation. Of course, within 
this region, the presence of the ghost is expected to signal instabilities. In particular 
one is tempted to interpret the failure of the full Vainshtein's resummation found in 
[TO] as linked to the presence of this sick DOF, even if one should be carefull to draw 
conclusions from energy arguments in General Relativity. Note that the presence of 
this ghost-like DOF is a generic feature of massive gravity as shown in the work of 
Damour and Kogan [5j. 

We also compared massive gravity with other models, namely deconstructed 
theories of gravity, as well as the DGP model. In the latter case we argued that the 
Vainshtein resummation process was of different nature, not involving a ghost degree 
of freedom. There are other variations on the DGP model which could free of ghosts 
too. For instance [3U] describes a model in which the five (or higher dimensional) 
Einstein-Hilbert action has a profile in the extra dimensions (or alternatively, has 
a varying Planck mass in the extra space). In this case, the absence of ghosts or 
strong coupling naturally introduces an incurable vDVZ discontinuity in the theory. 
We also note that the setup of "soft massive gravity" (31] uses higher dimensional 
generalizations of the DGP model where the usual induced Einstein-Hilbert action 
is supplied with additional UV operators, and [3T| argues that it does not suffer 
from any strong coupling issues. One interesting other way to avoid problems of 
massive gravity might be to break Lorentz invariance [321 HHH IHH IHHj • 
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Appendix: DOF reduction in the general case 

We start from the Lagrangian (JBJ) 

£ = ~9^<f> + ^(W) 3 - D0OW) 2 ), (64) 

and want to analyze the degrees of freedom it propagates. There are several ways 
to proceed. There is the "constructive approach" (see e.g. |2Sj, who use similar 
mechanism to reduce higher derivative gravity theories), or the formal Ostogradski 
method for higher derivative scalar field theories, as in e.g. [2U|. Both methods yield 
in fact the same result, and we will only discuss here the first one which is in direct 
correspondence with what we did in section 121 
One can rewrite the Lagrangian as 

C = ~d li </>& i <f> + F{d li d u <f>), (65) 

where F(X^) = js{X^ - 1^1^). Next we consider 

1 dF 

£ = "o^W + -tt^(cW - X^) + F ( x ^)- (66) 

The equations of motion for and <p obtained from this Lagrangian read 

dF 

= + - (67) 

= f F [d a dpcf>-X a p]. (68) 

Those equations reduce to those deduced from the Lagrangian (J54*)l if the 16 x 16 
matrix M^ v )^p) given by 

d 2 F 

M M , M = , (69) 

is invertible. In our case this matrix is given by 

<r 7 (3^7]^ X - rj^X^ - r]^X^ - Xrf a r) vP ) , (70) 



dX^ u dX af3 A 5 
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where X stands for X iliV if' v . The determinant of the above defined matrix reads 

det.M = -2 18 X 14 A- 80 (AX 2 + X aP X aP ) , (71) 
and does not vanish in general. Now we define 

*- = ^ < 72 > 

which we invert to get X^(tt) as a function of tt^ u . This inversion is guaranteed to 
exist again (al least locally) when M. can be inverted, but the explicit form of the 
inverse function may not be easy to obtain. Also one may have to divide the theory 
into different branches if the function is not globally well defined, in a similar way 
as what was done in section |21 In any case, our new Lagrangian is now 

£ = ~\d^<P + 7^(3^ - X^(vr)) + G(tv), (73) 

where G(u pa ) = F(X pu (n piJ )). The function G might now be regarded as a potential 
for the field 7i pa (together with n^X^n)). The investigation of the nature of 
the degree (s) of freedom in ir^ u in our case simplifies drastically because (unlike 
generic higher derivative gravity for example) we have a special mixing with a scalar: 
the coupling tt^ v 'd p d p <j) guarantees that there is only one component of ix^ v that 
propagates due to mixing with <f>. Indeed, decomposing the symmetric tensor into 
its orthogonal components 

^ = h& + + (iT - ^-)a + ^pA, (74) 

where and A v t are respectively transverse-traceless and transverse, we see that 
only A gets a kinetic term by mixing with <fi. The above form of the Lagrangian 
ensures that A is always a ghost. Indeed the Lagrangian is easily diagonalized by 
the transformation (f) — > ip = <fi — A. The final Lagrangian for the two DOF is: 

£ = ~\d^d^ + ^A^A - V(n pa (\)), (75) 

with V{jtpa) = ix^ X pv (y; pa ) — G(ir pa ). Although the above procedure is a bit formal, 
we did exactly this for the interaction (□</>) 3 , where indeed we encountered two 
branches and found a potential for the ghost A by inverting the scalar analogue of 
Eqn.(|72Jl. The exact form of the potential is now harder to find, having to solve 
equation (JT2*j) which reads 



' (3r]^X 2 - T]^X aP X al3 - 2XX^) (76) 



A 5 



for X^y as a function of 7r pcr . The least one can say by dimensional analysis is that 
the potential will be again IR relevant, as opposed to the UV relevant operator that 
we started from in (|64j) . 



19 



References 

[1] M. Fierz, Helv. Phys. Acta 12 (1939) 3; M. Fierz and W. Pauli, Proc. Roy. Soc. 
Lond. A 173, 211 (1939). 

[2] H. van Dam and M. J. G. Veltman, Nucl. Phys. B 22, 397 (1970). V. I. Za- 
kharov, JETP Lett. 12, 312 (1970) . Y. Iwasaki, Phys. Rev. D 2 (1970) 2255. 

[3] D. G. Boulware and S. Deser, Phys. Rev. D 6, 3368 (1972). 

[4] C. J. Isham, A. Salam and J. Strathdee, Phys. Rev. D 3 (1971) 867. 

[5] T. Damour and I. I. Kogan, Phys. Rev. D 66 (2002) 104024 
|arXiv:hep-th/0206042| . 



[6] G. Gabadadze and A. Gruzinov, arXiv:hep-th/0312074 



[7] C. Deffayet, G. R. Dvali, G. Gabadadze and A. I. Vainshtein, Phys. Rev. D 65, 
044026 (2002) |arXiv:hep-th/0106001| . 

[8] N. Arkani-Hamed, H. G eorgi and M. D. Schwartz, Annals Phys. 305, 96 (2003) 
|arXiv:hep-th/0210184| . 

[9] A. I. Vainshtein, Phys. Lett. B 39 (1972) 393. 

[10] T. Damour, I. I. Kogan and A. Papazoglou, Phys. Rev. D 67 (2003) 064009 
|arXiv:hep-th/0212155| . 

[11] R. Gre gory, V. A. Rubako v and S. M. Sibiryakov, Phys. Rev. Lett. 84, 5928 
(2000) |arXiv:hep-th/0002072| . 

[12] G. R. Dvali, G. Gabad adze and M. Porrati, Phys. Lett. B 484, 129 (2000) 
|arXiv:hep-th/0003054| . 

[13] G. R. Dvali, G. Gabad adze and M. Porrati, Phys. Lett. B 485, 208 (2000) 
|arXiv:hep-th/0005016| . 



[14] C. Deffayet, Phys. Lett. B 502, 199 (2001) |a rXiv:hep-th /0010186|. C. Def- 
fayet, G. R. Dvali and G. Gabadadze, Phys. Rev. D 65, 044023 (2002) 
|arXiv:astro-ph/0105068| ; 

[15] M. A. Luty, M. Porrati and R. Rattazzi, JHEP 0309, 029 (2003) 



|arXiv:hep-th /0303116|. 
[16] V. A. Rubakov, |a7Xiv:hep-th/0303125| 
[17] A. Nicolis and R. Rattazzi, |arXiv:hep-th/0404159| 
[18] G. Dvali, |arXiv:hep-th/0402l"30| G. Gabadadze, |arXiv:hep-th/0403161| 



20 



[19] C. Deffayet and J. Mourad, Phys. Lett. B 589, 48 (2004) 
arXiv:hep-th/0311124|. C. Deffayet and J. Mourad, Class. Quant. Grav. 
21, 1833 (2004) |arXiv :hep-th/ 03 11125] . 



[20 
[21 

[22 
[23 

[24 

[ 25 
t 26 



[27; 
[28 



[29 
[30 

[31 

[32; 

[33 
[34 
[35; 



F. J. de Urries and J. Julve, J. Phys. A 31, 6949 (1998) |arXiv:hep-th/9802115 



Teyssandier, P., & Tourrenc, P. 1983, Journal of Mathematical Physics, 24, 
2793 D. Wands, Class. Quant. Grav. 11, 269 (1994) [arXiv:gr-qc/9307034|. 



A. V. Smilga, arXiv:hep-th/0407231 



A. Hind awi, B. A. Qv rut and D. Waldram, Phys. Rev. D 53, 5597 (1996) 
|arXiv:hep-th/9509147| . 

R. Arnowitt, S. Deser and C. W. Misner, arXiv:gr-qc/0405109, 

S. A. Hojman, K. Kuchar and C. Teitelboim, Annals Phys. 96, 88 (1976). 

N. Arkani-Hamed and M. D. Schwartz, Phys. Rev. D 69, 104001 (2004) 
|arXiv:hep-th/0302110| ; M. D. Schwartz, Phys. Rev. D 68, 024029 (2003) 
larXiv:hep-th/0303114|; T. Gregoire, M. D. Schwartz and Y. Shadmi, JHEP 
0407, 029 (2004) | a rXiv:hep-th/0403224| . 

N. Kaloper,|arXiv:hep-th/0501028| N. Kaloper, Phys. Rev. D 71, 086003 (2005) 
larXiv:hep-th/0502035| . " 



M. Porrati, Phys. Lett. B 534, 209 (2002) | arXiv:hep-t h/0203014|. A. Lue, 
Phys. Rev. D 66, 04 3509 (2002) jarXiv:hep-th/0111168| . X~Gruzinov, 
|arXiv:astro-ph/0112246l T. Tanaka, Phys. Rev. D 69 (2004) 024001 
|arXiv:gr-qc/0305031| . 

G. Gabadadze and A. Iglesias, arXiv:hep-th/0407049 

M. Kolanovic, M. Porrati and J. W. Rombouts, Phys. Rev. D 68. 064018 (2003) 
larXiv:hep-th/0304148|; M. Porrati a nd J. W. Rombouts, Phys. Rev. D 69, 
122003 (2004) | arXiv:hep-th/0401211| . 

G. Gabadadze and M. Shifman, Phys. Rev. D 69, 124032 (2004) 
larXiv:hep-th/0312289| . 

N. Arkani- Hamed, H. C. Cheng, M . A. Luty and S. Mukohyama, JHEP 0405, 
074 (2004) |arXi v:hep-th/0312099| . 

V. A. Rubakov, [arXiv:hep-th/0407104| 



S. L. Dubovsky, JHEP 0410, 076 (2004) arXiv:hep-th/0409124 



G. Gabadadze and L. Grisa, arXiv:hep-th/04 12332 



21 



